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We describe in detail a counterexample to the topological classification of free fermion systems. 
We deal with a one dimensional chain of Majorana fermions with an unusual T symmetry. The 
topological invariant for the free fermion classification lies in Z, but with the introduction of inter- 
actions the Z is broken to Zg. We illustrate this in the microscopic model of the Majorana chain 
by constructing an explicit path between two distinct phases whose topological invariants are equal 
modulo 8, along which the system remains gapped. The path goes through a strongly interacting 
region. We also find the field theory interpretation of this phenomenon. There is a second order 
phase transition between the two phases in the free theory which can be avoided by going through 
the strongly interacting region. We show that this transition is in the 2D Ising universality class, 
where a first order phase transition line, terminating at a second order transition, can be avoided 
by going through the analogue of a high temperature paramagnetic phase. In fact, we construct 
the full phase diagram of the system as a function of the thermal operator (i.e. the mass term that 
tunes between the two phases in the free theory) and two quartic operators, obtaining a first order 
Peierls transition region, a second order transition region, and a region with no transition. 



I. INTRODUCTION 

The discovery of the quantum spin Hall effect Jl|, 
01 and of the strong 3D topological insulator 
both of which are novel band insulators, has prompted 
renewed interest in the study of topological phases of 
free fermion systems. Indeed, a full classification of all 
possible topological phases in such systems has been put 
forward in where it is related to the enumeration of 
symmetry classes of matrices , and in [llj , which uses 
the mathematical machinery of fC-theory. This classifica- 
tion is rather successful, with physical representatives of 
the non-trivial topological classes listed for dimensions 1 , 
2, and 3, including the quantum spin Hall system HgTe 
and the 3D topological insulator BiSb. 

The big open question now is how the presence of in- 
teractions changes this classification. Specifically, it is 
possible that phases that were distinct in the free classi- 
fication can actually be adiabatically connected through 
a strongly interacting region. Now, for certain systems 
the topological invariants can be defined in terms of phys- 
ically measurable quantities, and hence are stable to in- 
teractions. This occurs for example in the integer quan- 
tum Hall effect, where the integer Chcrn number is pro- 
portional to the Hall conductivity, as well as in 2D chi- 
ral superconductors and 2D topological insulators and 
superconductors. Also, the Z2 classification of the ZD 
topological insulator reflects the presence or absence of 
a TT theta term in the effective action for the electromag- 
netic field, extending the definition of this invariant to 
include systems with interactions [l^ [l^l • 

In this paper we give an example where the free classi- 
fication breaks down. The system is 1 dimensional, with 
an unusual T symmetry: = 1 instead of = ( — 1)^. 
For a concrete model, we consider the Majorana chain 
and its variations, where T acts on odd sites by TcjT~^ — 
—Cj so that terms like iCjCk are only allowed between 
sites of different parity. In the free-fermion setting, this 



symmetry is described by one positive Clifford generator, 
hence p=— 1, q = p -\- 2 ~ 1, and for d = 1 we get a 
topological invariant k e 7ro(i?q_d) = 

We can get some intuition for this integer by think- 
ing about boundary states. We start by comparing it to 
the usual Z2 classification of ID systems without sym- 
metry. The Z2 classification of systems without symme- 
try is reflected in the fact that for a pair of Majorana 
chains c", a = 1,2, we can gap out the dangling Ma- 
jorana operators c" and at the ends of the chain by 
introducing the terms iclcf and iclfcj.^ (the i is necessary 
to make these terms Hermitian). However, these terms 
are not T-invariant: T{iCjCj)T^^ — —iCjcj. Thus, in 

the T-symmetric case, we cannot gap out the dangling 
Majorana operators with quadratic interactions, for any 
number of chains - this is the origin of the Z invariant, 
which just counts the number of boundary states in this 
setup. However, it turns out that we can use non-trivial 
quartic interactions to gap out the dangling Majorana 
modes for the case of 8 Majorana chains - this is what 
we will focus on in this paper. 

Thus we study the setting of 8 parallel Majorana 
chains, which has a phase transition characterized by 
k = 8. We will see how the two phases separated by this 
transition are actually adiabatically connected through 
an interacting phase. This means that the two phases 
arc actually the same, and that the Z topological in- 
variant is actually broken down to Zg. Below we will 
demonstrate this fact by constructing an explicit path in 
Hamiltonian space connecting the two phases of 8 paral- 
lel Majorana chains. Adiabatic transformation along this 
path connects the two phases through a strongly interact- 
ing, but everywhere gapped, region. We do the analysis 
first for the microscopic model in section |TT1 where we 
construct a quartic T-invariant interaction that gaps out 
the 8 boundary Majoranas, and then for the continuum 
theory in section [1111 
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II. THE MICROSCOPIC MODEL 

We consider 8 parallel Majorana chains c", where a = 
1 , . . . , 8 is the chain index. The T symmetry still acts 
by Tcj'T~^ = —Sj. The Hamiltonian is the sum of the 
Hamiltonians for the individual chains: 



To motivate the construction, we must first delve a bit 
into some representation theory. In general, we can con- 
sider 2n Majorana's ci, . . . , C2n forming a 2" dimensional 
Hilbert space Ti (so in our case of interest n = 4). Ji 
is a representation of so{2n) as follows: given a skew- 
symmetric matrix A G so{2n), we let 
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^ u 



n-1 N 



(1) 



(2) 



(3) 
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The i in front is to make the matrix Hcrmitian, in order 
to obtain a unitary representation. It is easy to check 
that [—ip{A),—ip{B)\ = —ip{[A,B]), so that p defines a 
map of Lie algebras. The induced action on the Cj, 



We construct a path from a representative Hamiltonian 

of the u < V phase to one of the u > v phase. While in 
principle we could start with any representative Hamilto- 
nians for the two phases, it will be especially convenient 
to choose so-called fully dimerized ones, i.e. u = l,v = 
and u = 0,v = 1. This choice turns off the odd (or even) 
bond couplings and thus breaks down the chains into easy 
to analyze independent finite size systems. For example, 
for u = 1, t; = 0, the finite dimensional sub-systems con- 
sist of the Majoranas {cli_i,cli, . . . ,C2i_i,C2i}, and for 
u = 0,v = 1 they consist of {^f' ^2;+!, • • • , c|;, c^;_^ J. 
Both are 256 dimensional and we will generically denote 
their Hilbert space Hq. 

The key idea is now to work with these finite dimen- 
sional systems, in order to make a fully analytic treat- 
ment possible. Indeed, to connect the two phases, we 
start with one fully dimerized Hamiltonian, say u = 
l,v = 0, and turn on an interaction W which contains 
quartic terms, but only ones that are products of Majo- 
ranas with the same site index (but different chain in- 
dices), i.e. terms of the form c'j'^ c"^^ c"^^ c^"^ . The virtue 
of such an interaction is precisely that it does not couple 
the finite dimensional sub-systems; also, as required, it is 
T- invariant. We show that as we ramp up our specially 
constructed W, we can turn off the kinetic terms entirely 
(i.e. turn off u so that both u,v = 0) while maintaining 
a gap in Tio (and hence in the full system). We then 
reverse the procedure, turning on the opposite kinetic 
term v, and turning off W to get the opposite dimerized 
Hamiltonian u = 0,v = 1. 



Q i[p{A),ci] 



(4) 



A. Construction of the interaction term W 

We have described the path in Hamiltonian space that 
will connect the two phases; all that is left now is to 
construct the interaction term W in such a way that as 
W is turned on and the kinetic term turned off, the gap 
in Ho is maintained. We have thus reduced the problem 
to a completely finite dimensional one. 

W must couple different chains but only at the same 
site; thus it is made up of 8 Majorana fermions Ci, . . . , cg. 



is just the standard cj — J2j ^ij'^j- 

Before considering so(8), we warm up by studying 
so(4). Note that so(4) so(3) ® so(3). Under p, the 
generators of the two so(3)'s are 
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-(CiC2 - C3C4), -^{ClCA - C2C3), 2 (^1^3 + C2C4) ) (5) 



and 



1 z z 

-(C1C2 -h C3C4), -(C1C4 -I-C2C3), -(-C1C3 -I-C2C4) 

■ (6) 

The condition for a state \tp) G H to be annihilated by all 
generators of the first so(3) is equivalent to C1C2C3C4IV') = 
— and there are two states that satisfy this condi- 
tion, forming a spin- 1/2 representation of the other so(3). 
Similarly, C1C2C3C4IV') ~ IV') for lip) annihilated by the 
second so(3). Indeed, under so{3) (B so{3) , H decomposes 
as: 



W = (0,-) 



(7) 



From this analysis, we can already see that the above 
strategy for constructing the path through the space of 
gapped Hamiltonians connecting the two phases would 
fail for the case of 4 Majorana chains. This is because for 
4 Majorana chains the only possible quartic interaction 
W is proportional to C1C2C3C4, and regardless of its sign 
leaves a doubly degenerate ground state for each group of 
4 corresponding sites, leading to a gapless system at the 
midpoint of the path when the kinetic terms are turned 
off. Indeed, a field theory analysis shows that the corre- 
sponding marginal quartic interaction in the contimmm 
limit is a sum of two quartic terms, one for each so(3), 
but with opposite signs. One of these is irrelevant and 
the other relevant, leading to one of the so(3) sectors be- 
ing gapped out but the other remaining, producing the 
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(gapless) antiferromagnetic Heisenberg model. We will 
have more to say about the field theory analysis in the 
next section. 

The analysis of 4 chains does however lead to a natural 
way to gap out the 8 chain system. Indeed, we can split 
the 8 chains up into 2 groups of 4, turn on a quartic inter- 
action in each group of 4 to create low energy spin 1 /2's as 
in the above paragraph, and then couple these spin 1/2's 
into a non-degenerate singlet via an anti-ferromagnetic 
interaction. We will now flesh out this intuition and give 
an explicit construction of W. 

To facilitate the analysis, we make use of the so(8) 
symmetry. Per the construction above, we have 8 Ma- 
joranas forming a 16 dimensional representation so(8). 
This representation is actually a direct sum of two spinor 
representations, 8+ and 8_, of so(8), both of which are 
equivalent to the 8 vector representation in the sense that 
there is a group of so-called triality automorphisms of 
so(8) 14j which interchange the 8+, 8_, and 8. Now, 
while it turns out that a W that is fully so(8) symmetric 
doesn't work, we can find one that is symmetric under 
the so(7) C so(8) that fixes a particular element of say 
8+. This is the triality conjugate of one of the familiar 
so(7) subalgebras. 

The goal for the remainder of this sub-section is to 
show that the following explicit expression for W 

W = C1C2C3C4 -I- C5C6C7C8 -I- C1C2C5C6 

-I- C3C4C7C8 - C2C3C6C7 - C1C4C5C8 

+ C1C3C5C7 -I- C3C4C5C6 -I- C1C2C7C8 

- C2C3C5C8 - C1C4C6C7 + C2C4C6C8 

- C1C3C6C8 - C2C4C5C7 (8) 

is so(7) invariant. The proof that turning on W and 
turning off the kinetic terms leaves H gapped will be left 
for the next sub-section. 

To show that W is so(7)-invariant, we first have to 
specify the so(7) we are talking about. To do this, we first 
combine the Majoranas into regular fermions as follows: 



C2j-1 
C2j 



-i{aj - a]) 



(9) 
(10) 



We thus have 4 regular fermions. Let |0) be the state 
where all of them have occupation number 0, and let 
= -^(|0) — aJajflgaJlO)). 1-0) is a spinor in 8+, and 
we claim that it is the unique ground state of W , and 
that W is invariant under the so(7) subgroup of so(8) 
that leaves fixed. 

Let us start by proving is the unique ground state 
oiW . To do this, we first consider 



Wi = C1C2C3C4 -t- C5C6C7C8 + C1C2C5C6 -t- C1C3C5C7 (11) 

Note that each of the 4 terms in Wi has eigenvalues ±1. 
Re- writing the terms in Wi in terms of creation and an- 
nihilation operators, it is easy to see that is an eigen- 
state of eigenvalue —1 for all of them. It is also easy to 



see explicitly that jV') is the only state with this prop- 
erty - assuming eigenvalue (—1) for the first three terms 
immediately gives a state that is a linear combination of 
|0) and ajalagaJlO), and the correct linear combination 
yielding \'ip) is fixed by the last term. 

The key point now is that each term in W can be 
written as a product of terms in Wi - this immediately 
shows that is an eigenstate of each term of W , with 
eigenvalue ±1. In fact, it can be explicitly checked that 
all the eigenvalues are ~1, so that \'ip) is a ground state 
of W . Uniqueness follows from the fact that it is already 
a unique ground state of Wx . 

Now that we have proved that is the unique ground 
state of ly, we can finally show that W is invariant un- 
der the so(7) that leaves |'0) fixed. First, we identify the 
generators of so(7). Note that there are 28 linearly in- 
dependent bilinears CiCj. Take any one of the 14 terms 
in W, say CiCjCkCi, transposing a pair of c's if necessary 
to make the sign positive. We then claim that the set of 
all bilinears |(cjCj - CfeQ), §(ciQ - CjCk), ^{ciCk + CjCi), 
as CiCjCkCi ranges over the 14 terms in W, spans so(7). 
To see this, first notice that all these bilinears annihilate 
lip) - this follows from the fact that |0) is an eigenvec- 
tor of eigenvalue —1 for the corresponding terms in W. 
Thus the bilinears span a subset of so(7). To see that 
they actually span all of so(7) we use a dimension argu- 
ment - by brute force we compute the rank of the relevant 
28 by 28 matrix and find it equal to 21, the dimension 
of so(7). With this description of so(7), checking the 
invariance of W amounts to computing its commutator 
with all the generators of so(7). This is actually rather 
tractable, because for any given generator constructed in 
this paragraph, its commutator with all but 4 terms of 
W is trivially 0, and can be easily done via computer. 



B. Adiabatic continuation using W 

We now explicitly verify that we can connect the two 
phases adiabatically. As explained, this amounts to 
showing that we can turn on a quartic interaction and 
turn off the kinetic terms completely while maintaining 
a gap in the 256 dimensional Hilbert space Hq. That is, 
we have the 16 Majorana fermions c^, i = 1, . . . , 8, and 
consider the following interactions: 



(12) 
(13) 

tT (14) 



where W is as above, and W is given by replacing all 
the Ci's with cj's. Thus T is the kinetic term and W the 
quartic potential. 

We have explicitly diagonalized this Haniiltonian nu- 
merically. The eigenvalues are plotted in fig. [U where 
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we take the path w = 1 — t,t G [0, 1]. As we can see, the 
system remains gapped throughout. 
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FIG. 1: Eigenvalues oi H = tT + {1 — t) Wtot as a function of 
t. The system remains gapped throughout the path. 

However, it is nice to also have a clear analytical ar- 
gument that the system remains gapped. To do this, let 
us first analyze the symmetries. The 256 dimensional 
Hilbert space Hq for these 16 Majorana fermions has an 
action of the Lie algebra so(8)®so(8), with one so(8) act- 
ing on the Ci's and the other on the c^'s. The potential 
term Wtot is invariant under an so(7) © so(7) subalgebra, 
whereas the kinetic term T is invariant under the diago- 
nal so(8). In fact, to construct an easy to analyze Hamil- 
tonian, it will be useful to add in a fully so(8) © so(8) 
invariant quartic term Vtot — V + V' . V and V' are quar- 
tic in the fermions like W and W' , and are proportional 
to the quadratic Casimirs of the corresponding so(8)'s 
(plus some constant term); we will set the constant of 
proportionality below. The Hamiltonian H thus is 

H = tT + w Wtot + Wtot (15) 

The plan for the remainder of this section is to use the 
symmetries to reduce the problem of finding the spec- 
trum of H to something more manageable. We first 
study the actions of T, Wtot, and Vtot separately, and 
then block decompose H, with the blocks corresponding 
to different representations of the diagonal so(7), which 
is a symmetry of all 3 terms, and hence H. Then we 
construct a path through {t,w,v) space which connects 
(1,0,0) and (-1,0,0). 

A convenient basis for understanding better the action 
of T is as follows. We define bi, 6| through 

c^ - ih + bl) (16) 
= -lib. -hi) (17) 



We note that each term icjC- = blbi + bib] is equal 
to ±1, depending on the occupation number of the b 
fermion on site i. Thus the possible eigenvalues of T are 
±8, ±6, ±4, ±2 and 0. Under the diagonal so(8), the cor- 
responding eigenspaces have dimension 8-choose-fc where 
k is the number of occupied sites, i.e.: 1, 8, 28, 56 and 70. 
It will be convenient to consider superpositions of the m 
and — m eigenvalues of T and deal with representations 
li,,8i,,28L,56i and 1_r, 8i{, 28_r, 56_r, where L and R are 
exchanged under the action of T. All of these are irre- 
ducible, whereas 70 = 35^* © 35", where 35^ and 35" are 
distinct 35 dimensional representations. 

Now that we understand the action of T, let us figure 
out the actions of Wtot and Vtot ■ As a representation of 
so(8) © so(8), we can think of Hq as a tensor product 
of two 16 dimensional Hilbert spaces, one for the c^'s, 
and one for the c^'s (in fact, it is a graded tensor prod- 
uct, but this distinction will not make a difference in 
our analysis). This decomposition is useful because it is 
preserved by Wtot- As representations of so(8) we also 
have 16 = 8_ © 8+ where 8_ and 8+ are the two dif- 
ferent chirality spinors. These are distinguished by the 
sign of the fermion parity operator, with 8+ having even 
fermionic parity, and 8_ having odd parity. Indeed, we 
have two fermion parity operators, (—1)'^^ = Hi^i and 
(— 1)^" = nj=i '^'f well as the total fermion parity op- 
erator (-1)^ = (-l)-^^+-^^. The potential terms Wtot 
and Vtot respect both parities, whereas the kinetic term 
respects only (—1)'^, and flips (—1)''^^ and (—1)''^^. 

We can now expand the tensor product: 

Ho = (8-©8+)(g)(8_©8+) = (18) 
= (8_ «) 8_) © (8_ (E) 8+) © (8+ (g) 8_) © (8+ (g) 8+) 

With this description of the Hilbert space, it is easy to 
figure out the actions of Wtot and Vtot- Indeed, Wtot 
is a sum of two terms, W and W', that act indepen- 
dently on the left and right hand factors in the tensor 
product. Both W and W' are quadratic in the Lie al- 
gebra generators and commute with so(7) (which means 
they are proportional to the quadratic Casimir of so(7), 
plus a possible constant). We note that under the so(7), 
8-^8 and 8+ — > 1 © 7. Using these facts and doing 
some computation, we find that W annihilates 8_, has 
eigenvalue —14 on the state G 8+ that is fixed by 
the so(7), and has eigenvalue 2 on the remaining vector 
multiplet 7 of so(7), and similarly for W'. The actions 
of V and V, which are so(8) invariant, are even simpler: 
they assign a different energy to 8+ and 8_. We choose 
the coefficients so that this energy is equal to for 8+ 
and 1 for 8_. 

We now want to relate the two descriptions ^TE\i and 
(UHl) of Ho to find the action of H = tT + w Wtot + 
V Vtot ■ As a first approximation to a connecting path, we 
attempt to connect the phases in a purely so(8) invariant 
way. That is, we connect the points (1, 0, 0) and (—1,0, 0) 
in the {t, w, v) space of Hamiltonians (I15p by varying only 
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t and V. Specifically, the path is 

(t, w,v) = (cos 61, 0, sin 61) , e [0, tt] (19) 

The computation of the spectrum along this path re- 
duces, by so(8) symmetry, to diagonalizing the blocks 
corresponding to the 1, 8, 28, 56, and 35 representations. 

Now, we have the following decomposition of the tensor 
factors in (|18p as representations of the diagonal so(8): 

8_ ®8_ = 10 28 ©35^ (20) 

8_ ® 8+ = 8 56 (21) 

8+0 8_ = 8 56 (22) 

8+®8+ = 10 28 35" (23) 

Let us figure out the action of T. By looking at the 
fermion number, we see that T exchanges ([20| and ([23]) . 
and (HD) and (j^ . The eigenvalues of T on the pairs 
of 1, 8, 28, 56, and 35 representations are ±8, ±6, ±4, ±2, 
and 0, and v V just assigns an energy v to ([2T|) and (p2)l . 
and an energy 21; to ([20|) . So is represented by the fol- 
lowing 2 by 2 matrices on these pairs of representations: 



From the discussion above, we know that the action of 
Wtot on this decomposition is as follows: Wtot has eigen- 
value ~28w on 1 (g) 1, — 12w on 1 (g) 7 7 (X) 1, and Aw on 
7 7. To relate this decomposition to we first use 
the fact that under so(7), 28 7021, and 35 107027 
dl. Now, 7®7= 1021027. Let H = V + tT + w Wtot- 
We immediately have the eigenvalues 4w on the 27 and 
Aw — 8i^ on 21. Also, the two 7's have the eigenvalues 
— 12w and —12w — 8t^. All that is left is to compute the 
eigenvalues of H on the two I's, and show that one of 
them is lower than the other and any of the above. 

To do this, let us denote the 1 in 8+ = 1 7 by \ip), 
and let \ip'^) = \ip) ^ \tp)- Then the so(8) invariant 1 in 
8+ 8_ can be written as 



|U(8)) = -^ (29) 

Let P = |lso(8))(lso(8)l- The effective Hamiltonian on 
the two dimensional space of I's is then 



l-( ^ (24) 
\8 cose 2 sin6i J ^ ' 

S ) (25) 

\^ 6 cos smd J ^ ' 

28: , (26) 
\ 4 cos 2 sin& I ^ ' 

_„ f sm6 2 cos6' \ 
56 : „ . ■ Q ] 27) 
\^ 2 cos smO J ^ ' 

with the two 35 's splitting up into the 1 by 1 matrices 
2 sin^ and respectively. Diagonalizing these, we see 
that the lowest eigenvalue occurs in the 1 sector; it is 
sin 9 - ^8 cos2 6 + sin^ 9. 

The problem with this path is that there is a large 
degeneracy at 9 = tt/2, where the lowest eigenvalues for 
the 1, 28, and the 35 representations are 0. To break this 
degeneracy, we will turn on the so(7)-symmetric Wtot 
term. We can work locally around 9 = tt/2, where (with 
a suitable rescaling) we can set w = sin w 1 , and keep 
t = cos 9 infinitesimal. With such infinitesimal t, the 
space of low energy states is well approximated for our 
purposes by 8+ (g) 8+ = 1 28 35, see ([23]). To second 
order in t, the eigenvalues oiV + tT are —32^'^ on the 1, 
—St"^ on the 28, and on 35. To understand the action 
of Wtot on 8+ (g) 8+, we note 

8+(g)8+ = (10 7)® (10 7) 

= 1 «) 1 (1 ® 7 7 (g) 1) (7 ® 7). (28) 



H,s = -32fP-28w\^p^){ij^\+'iw (1 - \i>^){'4'^\) , (30) 

which is represented in matrix form by 

/-28u;-4<2 -4V7t2 \ 

-477^2 Aw ~ 28e ) ■ ^ ' 

We graph the eigenvalues of this matrix, as well as those 
for the 7, 21, and 27 representations, in figure [2l 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 2: Eigenvalues of H + tT + w Wtot along the path 
(tjw) — (0.1 cos ^, 0.1 sin ^), S £ [0,7r]. The degeneracy is 
broken and the system remains gapped. 
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III. CONTINUUM ANALYSIS 

We now turn to a field theory analysis of the above 
phenomenon. In the free system, the transition is simply 
that of 8 parallel decoupled Majorana chains, so its field 
theory description is 

iJfrco = Hq + Hjn (32) 

= 2 ^ / ^^^^^^^ ~ ^j^^h) dx (33) 
+ ^ XI / ^VjVjdx. (34) 

Here the mass m is a free parameter, and the phase tran- 
sition occurs at m = 0, where the system is described 
by the conformal field theory (CFT) of 8 free Majorana 
fermions, or, equivalently, the 6*0(8)1 WZW model. We 
have shown in the previous section how to smoothly con- 
nect the phase with positive m to one with negative m 
through a strongly interacting region in a microscopic 
lattice model. In this section, we construct a field theory 
analogue of that path. 

To study the field theory of the transition, we will need 
to understand the structure of the critical point in more 
depth. As we said, this is just the CFT of 8 free Ma- 
jorana fermions, which, in addition to the fermions also 
contains spin fields, and these will be useful in our con- 
structions. For the case of one free fermion, there are the 
conformal weight (1/16, 1/16) order and disorder opera- 
tors a{z,z) and fj,{z,z). For the case of 8 free fermions, 
we can form the 2^ = 256 possible products of ai and 
{i = 1, ... ,8). These have conformal weights (1/2, 1/2) 
and are precisely (linear combinations of) the 16^ pri- 
mary fields -ipi-ipj, xSj, i^iXj, XiXj, hj = ^,---,8, where 
■07 X? "0: X the two chiral left and right moving spinor 
representations of S0{8) - indeed, this is a manifestation 
of the equivalence between the WZW 50(8)1 and free 
fermion models p^ . 

The 50(8)1 WZW model has a triality symmetry 
which permutes the rj, ip, and x fields. The fact that 
the Hamiltonian (j32[) is invariant under triality, as well 
as other facts we will need, can be seen from bosonizing 
the system. We pair up the 8 Majorana fermions rjj into 
4 Dirac fermions and bosonize those: 



r]2j-i±ir]2j = ex.p{±i(f)j) (35) 

fj2j-i±i'i]2j = exp{±i4)j) (36) 

The 16 (1/2,1/2) fields 

exp^(±0i±02±03±04) (37) 



are then the spinors tJj and x, with chirality distinguished 
by the parity of the number of minus signs in ([37]). By 



doing a change of coordinates in space, we can go to 
a basis in which the spinors look like the vectors, etc. - 
indeed, we can implement the entire triality symmetry 
group by appropriate rotations in space. 

Now that we have an understanding of the critical 
point, we look at perturbing it with the quartic inter- 
action W ^ from the previous section. Let us first de- 
termine the continuum limit of this interaction. From 
the action c* — + (— l)^c* of the T symmetry we see that 
at low energies 

c^fe = 7y,(2fc)+j?,(2fc) (38) 
C2k^i - mi^k - 1) - f,,{2k - 1) (39) 

with T just swapping rji and fji. Now, W is a sum of 
14 quartic terms; to start off we analyze just one of 
these, and we will take the sum over two adjacent sites: 

^2k-1^2k-1^2k~1^2k~l + ^2k ^2k ^2k ^2k- We have 

C2fe-lC2fe-lC2fc-lC2fe-l + C2fe C^^. C^^. C^k = (40) 
4 4 

X{{^l^+f]i) + \{{r1,~f]i) (41) 

i=l 1=1 

Expanded out, this gives 8 terms, 6 of which are nonchi- 
ral, with the other 2 being products of all ?7's or all fy's. 
Instead of trying to argue that these last 2 are irrelevant, 
we will simply ignore them - after all, we have the free- 
dom to choose whatever continuum Hamiltonian we wish. 
The remaining 6 terms can be reorganized suggestively 
as follows: 



2 [vim + mvi) (mm + mm) (42) 

{viV2 - mv4:) {mm ~ mf]4) (43) 

+^ {viV4: + mm) {mm + '72%) (44) 

{mm -mm) {mm -mm) (45) 

+^ {mm - mm) {mm - mm) (46) 

{mm + mm) {mm + mm) (47) 

This can be compactly rewritten as 

3 3 

^E^^-^a-^E^a^a (48) 
a—1 a—1 



with Ja and being the generators of the two SU{2) 
groups in the 50(4) of the 4 chains. 

Now, when we take the sum over all 14 such quar- 
tic terms in W, we get a more complicated expression. 
However, we know that it is a sum of terms of the form 
mVjViVjj i-S- a sum of products of right moving and left 



moving so(8) currents. We also know that it is so(7) 
invariant, where the so{7) leaves a spinor -^g fixed. In 
addition, we know that it is T invariant. Using the fact 
that the so(8) currents can be equally well expressed as 
bilinears of the spinors ipj (this follows from bosonization 
([35l [37|) ). we see then that the only allowed interactions 
are of the form 



int 



- B 



M& (49) 



This follows from the fact that, with the T symmetry, the 
allowed interactions are quadratic forms on so(8). As 
a representation of so(7), so(8) = so(7) 7, and both 
of these irreducible factors have exactly one invariant 
quadratic form. 

For the remainder of the analysis, it will be convenient 
to use triality to re-express the free Hamiltonian ((32|) in 
terms of the ij:, i.e. treat the as the fundamental fields. 
Without a mass term, it is: 
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E 



{ijjjd'iljj — ipjdil^j) dx 



(50) 



We now view both (|49)) and the mass term m^,^^ rjjfjj 
as interactions, and construct a 3 dimensional phase di- 
agram in A, B, and m (the parameters A, B, and m will 
have non-trivial renormalization group flows, and for the 
purposes of the phase diagram are all defined at some 
fixed energy scale). To do this, we first need to express 
the mass term in terms of the free fermions ipj. If we 
look at the ipj as fundamental vectors, the r]j and fjj are 
then spinors and can be expressed as linear combinations 
of products of order and disorder operators. Comparing 
the bosonizcd form of the order and disorder operator 
bilinears [13, Ei 



' / n " n 



(54) 



The goal for the remainder of this section will be to show 
that we can smoothly connect the region to > to the 
region to < by turning on A and B. To accomplish this, 
we first work out some features of the phase diagram - 
the full phase diagram is shown in figure ([5]) . We can get 
some information about the phase diagram by studying 
specific special points. Indeed, the fine B = 2A, {A > 0) 
gives the marginally relevant interaction 



= -A 



(56) 



This is the Gross-Neveu model, which is known to be 
integrable and gapped, and in which (riifjj) = M6ij con- 
denses, spontaneously breaking the Z2 chiral symmetry 
'^3^3 ~^ ~'^3^3 (indeed, the spectrum of the model and 
the action of triality on the different types of excitations 
is well understood - see [l3|)- The two vacua have oppo- 
site signs of M, and are adiabatically connected to the 
+171 and —TO gapped phases respectively. Thus, to con- 
nect the +m and —to phases, all we have to do is connect 
these two vacua, which we can think of as corresponding 
to infinitesimal points to = ±e. 

The line B = —2A, {A > 0) gives another Gross-Neveu 
point, related to the previous one by the chiral transfor- 
mation V-'s ^ ""08 J "08 ~^ "08, and the lines B = 2 A 
and B — —2A for A < give rise to marginally irrele- 
vant interactions that flow to the gapless SO {8)1 theory. 
Also, the line B = gives rise to the so(7) symmetric 
Gross-Neveu model 



CT2fc-i cr2k = sm 



A*2fe-i A*2fc = cos ■ 



0fc - 4'k 

2 _ 

4'k - (t>k 



(51) 



(57) 



and the bosonization of the mass term via ([551 \T7} we 
can show that (up to a constant factor) 



■ n + n 

3=i 3=1 



(52) 



Note that here the Cj and fij are order and disorder op- 
erators for the tjjj fermions. 
The total Hamiltonian is thus 



^ = ^ e/ (^^^^^- 



Ipj dip j) dx 



(53) 



This line provides the most important insight into the 
problem. It has a gapless Major ana mode, ipa and we can 
explore its neighborhood by turning on an inflnitesimal 
B: 



Hint = I E "^^^1 - ^ E ^J^J- I ^8^8 (58) 

y3 = l J \3 = 1 

with B <^ A. We can solve this Hamiltonian by treating 
the B term as a perturbation. Thus we have an expecta- 
tion value V'i'^j = -^'j ^'^^ the effective low energy 
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theory is simply that of the free fermion -08 with mass 
BM': 



(59) 



- m {{a)' as - (/i>^ Ms) - BM' 08^8 (60) 

There are two possibihties: if the condensate M' = 
V'jV'i > then the aj,j ~ 1, . . . , 7 gain expecta- 
tion values and m couples to the Ising order parameter 
tTgi whereas if M' < 0, then Mj, j = 1, . . . , 7 gain expecta- 
tion values and /zg plays the role of the order parameter. 
Either way, the low- and high-temperature phases of the 
effective Ising model are realized at i? > and B < Q 
respectively. 



m 



71 

/ m= e , B=2A 



^ B 



\ m=- e , B=2A 



Ising transition 
region 



FIG. 3: Smoothly connecting the m > and m < phases 
through an Ising transition region 

We now have enough information to connect the two 
phases - see figure ([3|). First of all, the two phases are 
connected to the m — zLe points on the B — 2A, A > 
line. From the RG diagram (see discussion below), un- 
less something unexpected happens at strong coupling 
we see that these points are adiabatically connected to 
the points < B ^ A, m = ±e. But here we can 
work with the effective Hamiltonian If we go from 

(_B, +e) to (-B, — e) by going around the origin in the 
{B,m) plane, through the B < region (so as to avoid 
the positive B axis), the quantity Z — 5X]j=i'?j'?j ~ 

— 11^=1 '^j + n^=i Mj (see ([52]) ) will vary smoothly. This 
is the desired path that connects the two vacua. Thus, 
because B controls the thermal operator and m the order 
(or disorder) field, the critical theory turns out to be in 
the universality class of the 2D Ising transition. 

For completeness we now work out the entire renor- 
malization group flow diagram. The two operators 




(61) 
(62) 



form a closed set under RG flow, as they are the only 
so(7) invariant (1, 1) operators available. To track the 
flow of the coupling constants A and B we use Polyakov's 
formula [Ti 



dgn 



dln(A/fc) 



(63) 



where the C^^ are the coefficients of the three point func- 
tion of the corresponding operators. Thus we are reduced 
to computing 3 point functions of products of X and Y. 
This is trivial, since we are in a free fermion theory and 
can apply Wick's theorem. We immediately see that, 
since Y is the only operator containing , we must have 
an even number of Y^s for a nonzero coefficient. So the 
only nonzero coefficients are C-^^ , Cy^ , Cy 
The RG flow equations then read 



■■Y ' 



and CY ■ 



dA 



d\n{K/k) 

dB 
d\n{A/k) 



SA^ + TB^ 
UAB 



(64) 
(65) 



for some constants S,T,U > 0. Because B ~ 2A is a 
fixed line, U = S + AT. We do not calculate S and T; 
all positive S and T give qualitatively similar results for 
the RG flow diagram, which is reproduced in figure (j4]) 
below. 



B=-2A 



8=0 



B=2A 




FIG. 4: RG flow for B (a;-axis) and A (i/-axis). 

From the flow diagram fS]), we can construct an entire 
phase diagram for the transition from positive to negative 
m. Indeed, all points with B > that are above the line 
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B=-2A B=0 B=2A 




transition 
region 

FIG. 5: Phase diagram for the transition from negative to 
positive m, indicating dependence on B (a;-axis) and A (y- 



B = —2A flow to the same fixed point and correspond to 
a first order transition. As we saw, there is no transition 
in the B < region, and this extends to all points above 
the line B — 2A. Finally, the remaining points flow to 
the origin, which corresponds to free fermions and thus 
indicates a second order transition. The situation is sum- 
marized in figure 
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